Adiabatic transport of Bose-Einstein condensate in double-well trap 
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A complete irreversible adiabatic transport of Bose-Einstein condensate (BEC) in a double-well 
trap is investigated within the mean field approximation. The transfer is driven by time-dependent 
(Gaussian) coupling between the wells and their relative detuning. The protocol successfully works 
in a wide range of both repulsive and attractive BEC interaction. The nonlinear effects caused by 
the interaction can be turned from detrimental into favorable for the transport. The results are 
compared with familiar Landau-Zener scenarios using the constant coupling. It is shown that the 
pulsed Gaussian coupling provides a new transport regime where coupling edges are decisive and 
convenient switch of the transport is possible. 



I. INTRODUCTION 

Nowadays the trapped Bose-Einstein condensate 
(BEC) is widely recognized as a source of new fascinat- 
ing physics and remarkable cross-over with other areas 
0, 0, y, 0, H, H, 0| ■ Between diverse directions of this 
field, a large attention is paid to dynamics of weakly 
bound condensates (or multi-component BEC) with em- 
phasis to nonlinear effects caused by interaction between 
BEC atoms. The studies cover different aspects of bo- 
son Josephson junctions B, M , M , [ToL fill [Hj], popula- 
tion ^of topological states 113 . 03] > tunneling processes 
18 . Jl9l |20. [211 ] , t ransport of condensate 
"lEWlaliSlslSi and other topics. Both 




multi-level or multi-well systems are considered. In the 
first case, BEC is confined in a single-well trap and con- 
tains atoms in a few hyperfine levels, thus forming BEC 
components. The laser light can couple the components 
and initiate various regimes of the population transfer. In 
the second case, BEC atoms occupy the same hyperfine 
state but the trap is separated by laser-produced barriers 
into a series of weakly bound wells. The atoms can tun- 
nel through the barriers and exhibit similar effects as in 
a multi-level system. In this case BEC components are 
represented by populations of the wells. 

The present work is devoted to complete irreversible 
transport of BEC in a two-well trap. We will consider 
the process when BEC, being initially in one of the wells, 
is then completely transferred in a controllable way the 
target well and kept there. Such transport can be re- 
alized in multi-well traps [l(| and arrays of selectively 
addressable traps [3(J. Being useful for a general ma- 
nipulation of the condensate, the irreversible transport 
can also open intriguing prospects for producing new dy- 
namical regimes [20], investigation of topological states 
[HI EH ' generation of various geometric phases [2(| [3l| . 
The latter is especially important since geometric phases 
are considered as promi sing information carriers in quan- 
tum computing [32], |33j, |34| . 

In general the BEC transport can be designed by var- 
ious ways: Rabi switch time-dependent potential 




modulation [22|, Rosen-Zener 
batic population transfer [H, |24|, |25|, |26|, |27|, |28 
Adiabatic methods seem to be especially effective for the 
transport because of their robustness to modest varia- 
tions of the process parameters [H, [H, [37| ■ In particular, 
the Stimulated Rapid Adiabatic PassageVsTIRAP) [Hj 
is widely used for this aim @, [H HlUS HHH- This 
method promises robust and complete population trans- 
fer of the ideal condensate but suffers from nonlinearity 
effects caused by interaction between BEC atoms. Then 
a natural question arises, is it possible to turn the nonlin- 
earity from detrimental to favorable factor of adiabatic 
transport? 

In this respect the study of nonlinear Landau-Zener 
(LZ) tunneling of BEC in the systems like two-band ac- 
celerated optical lattices give a useful and promising mes- 
sage 0, [HI, [2(| ■ It was shown that in such systems the 
tunneling is asymmetric, i.e. it can be considerably en- 
hanced or suppressed by the nonlinearity, depending on 
its sign. This means that the LZ scheme [39[ can serve 
as a good basis for nonlinear transport of BEC. 

The present paper is devoted to realization of this idea 
for adiabatic transport of interacting BEC between two 
potential wells. In addition to the familiar LZ case with 
a constant coupling Q between the wells [39(, we will in- 
spect the nonlinear transport with time-dependent cou- 
pling fl(t) of a Gaussian form. Actually this a general- 
ization of both LZ (constant coupling and linear bias of 
diabatic energies) [391 ] and RZ (time-dependent coupling 
and constant diabatic energies) [40( methods. Such gen- 
eralization is partly motivated by Stark Chirped Rapid 
Adiabatic Passage (SCRAP) method [il|. Being based 
on LZ protocol with a pulsed laser coupling, SCRAP 
turned out to be quite effective in atomic physics. Note 
also that a constant coupling is somewhat artificial and 
its time-dependent version is more realistic. 

As shown below, the nonlinear transfers with constant 
and time-dependent couplings have much in common 
(e.g. asymmetric impact of nonlinearity, determined by 
the interaction sign) and both them can result in effec- 
tive adiabatic transport at strong nonlinearity. However 
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the time-dependent (pulsed) protocol is more flexible and 
has interesting peculiarities. For example, we get a new 
transport regime where the coupling edges play a deci- 
sive role. Besides, the pulsed coupling allows us to switch 
on/off the transport in much wider interval of the detun- 
ing rate a and control the switch value a s through the 
coupling width T. 

In the present paper these peculiarities are studied by 
using the stationary spectra, in particular the loops aris- 
ing there due to the nonlinearity. Although the transport 
problem is not stationary and we actually deal with time- 
dependent Gross-Pitaevskii equation [421 . such analysis 
is known to be very instructive [HI, Il6l.ll7l |20| . It gives a 
clear intuitive treatment of quite complicated processes 
and allows for getting simple analytical estimations. 

The paper is outlined as follows. In Sec. 2 the mean- 
field formalism is outlined. In Sec. 3 we present the 
transport protocol. In Sec. 4 results of the calculations 
are discussed. The summary is done in Sec. 5. 



II. MODEL 

A. General formalism 

The calculations have been performed in the mean- 
field approximation by using the nonlinear Schrodinger, 
or Gross-Pitaevskii equation [42j 



2m 



-V ext (r,t)+g \^(r,t)\ 2 ]^(r,t) (1) 



where the dot means the time derivative, ^(r, t) is the 
classical order parameter of the system, V ex t t) is the 
external trap potential involving both (generally time- 
dependent) confinement and coupling, go = Aita/m is 
the parameter of interaction between BEC atoms, a is 
the scattering length and m is the atomic mass. 

The order parameter of BEC in a double- well trap can 
be expanded as Q 



(2) 



= ViV(V>i(i)$i(r) + iMt)S 2 (0) 



where $fc(f) is the static ground state solution of ([I]) for 
the isolated k-th well (k — 1,2) and 



(3) 



is the corresponding amplitude expressed via the relative 
population N k (t) and phase <j) k {t). The total number of 
atoms N is fixed 



dr\V(r,t)\ 2 /N = JVi(t) + N 2 (t) = 1 



(4) 



Substituting @ to (JT|) and integrating out the spatial 
distributions we get Q 



#i = [Ei(t) 
#2 = [E 2 (t) 



LW|V>i| 2 ]V>i 
IW|V> 2 | 2 ]V> 2 



(5) 
(6) 



where 

is the coupling between BEC fractions 

Ek {t) = iy dr[^-\vn\ 2 + nv ext m k ] 

is the depth of the k-th well, and 



(7) 



(8) 



(9) 



labels the interaction between BEC atoms in the same 
well. Here we assume that U\ = U 2 = U '. The values 
fl(t), E k (t), and U have the dimension of frequency. 

In the present study we use a Gaussian coupling with 
a common amplitude K 

Q(t) = KQ(t), n(t)=«p{-fc^} (10) 

where t and T are centroid and width parameters. Then 
scaling ©-© by 1/(2K) gives 



#1 = [Bi(t)+A|V>i| 2 ]Vi-^(*)^2 

1^2 = [E 2 {t)+k\^ 2 \ 2 ]i>2-~&{t)^l 



where 



E k (t) = 



E k (t) 
2K ' 



A 



UN 
~2K 



(11) 
(12) 

(13) 



and the scaled time, 2Kt — ► t, is dimensionless. In (fTi"]) - 
(fT2"]) . the key parameter A is responsible for the interplay 
between the coupling and interaction. 

Substituting ([3]) into (fTTj) - (|12p leads to the system of 
equations for the populations N k (t) and phases <p k {t): 



N k = -Q(t)^/NjNk Bin(0j - fa) 



h = -[E k (t)+AN k ] + -Q(t) 




(14) 



with j ^ k. Considering N k and — <p k as conjugate vari- 
ables and using the linear canonical transformation (26j 



z = Ni-N 2 , 

.-if 



,Z = Nx + N 2 = l 



»2 



) ,e = -i 



i 



+ 4>2) 



(16) 
(17) 



one may extract the integral of motion Z and correspond- 
ing total phase <d from equations ([TJJ-tlS]) and thus re- 
duce the problem to a couple of equations for new un- 
knowns, population imbalance z and phase difference 9. 
These equations read 



-Q(t)\/l - z 2 sin 29 
I[A(i)+A Z + fi(i)-^=L= 



cos 29} 



(18) 
(19) 
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where 



A(i) = £?i(i) -E 2 {t) =at. 



(20) 



is the detuning (difference in well depths). In accor- 
dance with LZ practice we assume for A.(t) a linear time- 
dependence with the rate a. 

Equations of motion (fT8| - (fT9|) are invariant under 
transformations 



A-f-A, a^-a, 9^-9 + ^ (21) 



a->-a, z^-z, 9^-9. (22) 

The former relates the transport protocols for repulsive 
and attractive BEC in one direction while the latter con- 
nects the transport in opposite directions (with corre- 
sponding interchange of the initial conditions). 

Note that equations (fT8|) -(fT9 | allow a classical anal- 
ogy with z and 9 treated as conjugate variables. It is 
easy to verify that these equations can be recast into the 
canonical form 



dRci 
89 



9 



with the classical Hamiltonian 

H d = ^[A(t)z + ^z 2 - Cl{t)y/l - z 2 cos 26»] 
and chemical potential [43j 



(23) 



(24) 



\l = R c l + V ir 



i[A(t)z + Az 2 - fl{t)y/T 



z 2 cos 29} 
(25) 

where V int = Az 2 /4. 

It is easy to see that, up to notation, Eqs. (TT8|) . ([T9"]) 
and coincide with those in Q. However, unlike the 
previous studies of the oscillating BEC fluxes in traps 
with constant parameters @, [l2j], we will deal with ir- 
reversible BEC transport by monitoring time-dependent 
parameters E k {t) and O(i). 



B. Stationary states 

Since we are interested in adiabatic (very slow) trans- 
port, the analysis of stationary states can be useful. 
These states are defined by the condition 



z = 6 = 



(26) 



Despite the fact that our Hamiltonian and thus the vari- 
ables z and 9 actually depend on time, this dependence 
is assumed to be slow enough to make (J2HI) relevant. 
Under condition (f26|) equations (fT8| -(|19 p give 



(27) 
(28) 



7T 

= 2 U 



A(t)+z(A± 



VT^ 2 



) = o 



where n is an integer number. Equation (|28[) has "+" 
or "— " for even and odd n, respectively. The analysis of 
(|28f finds two real roots for |A| < and four real roots 
for |A| > fj, i.e. at high nonlincarity. 

Substitution of the stationary solutions of (|27[) - (j2"5|) 
into ([25]) gives the chemical potentials /J,-(t) and n+(t), 
representing the eigenvalues of the stationary states, i.e. 
(JHD casts into ip k (t) = ' N k {t)e~ ilMt . In the four-root 
case, the repulsive interaction (A > 0) leads to one so- 
lution for /z+(i) and three solutions for /z_(i), and vice 
versa for attractive interaction (A < 0). 

As shown below, three roots form a nonlinear loop in 
the stationary spectra. For |A| > fj, the time length of 
the loop reads [ijj 



tr 



(|A| 



2/3 



(fi(t c )) 



2/3\3/2 



(29) 



where f2(i c ) — > 1 for the constant coupling. It is seen 
that the length increases with magnitude of nonlinearity 
A and decreases with detuning rate a. The nonlinear 
structures (single loops, double loops, butterfly, ...) can 
appear in the stationary spectra for smaller nonlinearity 
|A| < £l as well, see discussion in Ref. [44l |. 



III. TRANSPORT PROTOCOL 

In the standard LZ protocol without nonlinearity 
(39| . the final probabilities of two competing processes, 
Landau-Zener tunneling between adiabatic levels (Plz) 
and adiabatic following (P), read 



P, 



LZ 



P = 1 - P, 



LZ 



(30) 



Both processes are controlled by a constant coupling f2 
and detuning rate a. A complete adiabatic transfer 
P = 1 takes place in the adiabatic limit a — > 0. As shown 
in [TU, [HI, [l3] , the nonlinearity caused by interaction be- 
tween BEC atoms essentially modifies the transfer (|30p . 
Namely, the diabatic LZ tunneling becomes possible even 
at a — > 0. Moreover, the impact of nonlinearity is asym- 
metric, i.e. it favors or suppresses the LZ tunneling and 
inversely affects the adiabatic following, depending on 
the interaction sign [lj|. This suggests a principle pos- 
sibility to use the LZ scheme for the complete adiabatic 
transport in the nonlinear case. 

In the present study, this point is scrutinized for a 
constant and time-dependent Gaussian coupling. The 
relevant protocols are illustrated in Fig. 1 where the 
detuning (f2T)]) is modeled by 



E^t) = ^at, E 2 {t) = —at 



(31) 



with a > and — t s < t < t s . We assume that at early 
time t < —t s the diabatic energies are tuned from the 
symmetry values Pi = E 2 = to the values Pi(— t s ) 
and E2(—t s ), then exhibit the linear evolution (13~1]) at 
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Constant coupling 

A=0 



Gaussian coupling 

A=0 




Constant coupling 

a) a = 



Gaussian coupling 
b) 



FIG. 1: a): Stationary adiabatic levels /x± (bold curves) 
calculated in the linear case (A=0) for the constant coupling 
Q=l (dotted line) and a=0.35. The diabatic levels _Ei ; 2 (strait 
dash curves) are given for the comparison. BEC is initially 
placed at the lower level as indicated by the bold dot. The 
arrows show the diabatic Plz and adiabatic P transfers with 
the probabilities (|30[) . b): The similar scheme for the time- 
dependent Gaussian coupling (I10|l . 



— t s < t < t s and finally come back to the zero symmetry 
values. In the experimental setup such manipulations can 
be produced by varying the well depths. In what follows, 
we will present the results only for — t s < t < t s . 

The conventional LZ scheme is illustrated in Fig. QJ,). 
Despite its artificial properties (constant coupling and in- 
finite diabatic energies at t — > ±00), it is known to give 
quite realistic transition probabilities (|30|) . The reason 
is that the LZ transfer actually occurs only during a fi- 
nite time of close approaching the diabatic energies near 
t ~ 0, when A(t) < f2. The lateral regions of early and 
late time, where the coupling impact is negligible, are 
irrelevant. Hence the LZ artifacts are not essential. 

However, if anyway LZ works only for a finite time, 
then it is natural to use a time-dependent coupling of 
a certain duration. In this connection, we propose an 
effective and flexible transport protocol where a time- 
dependent coupling fi(i) of a Gaussian form (fit))) is used, 
see Fig. QJ)). Then, in addition to the usual LZ control 
parameters, detuning rate a and coupling amplitude K, 
we get a new one, the Gaussian width T. As discussed 
below, this protocol allows a rapid and complete switch 
of adiabatic BEC transport in a wider range of a and pro- 
vides a new transport regime where not the center but 
edges of the coupling f2(£) are important. This makes 
the new protocol more effective and flexible. Hereafter 
it is referred as a Gaussian Landau-Zener (GLZ) proto- 
col. Obviously, the GLZ is also a generalization of the 
Rosen-Zener scheme [4fj| | which exploits a time-dependent 
coupling but constant diabatic energies. 

In the present calculations we use the Gaussian cou- 
pling (|10p with the centroid at i =0 and width param- 
eter F=5, apart from of Fig. 7 where various T are 
applied. The coupling full width at half maximum is 
FWHM=2.355 T = 11.78. The monitoring of the cou- 
pling Q(t) which is the penetrability of the barrier be- 
tween the wells, can be produced in experiment by vari- 
ation of the spacing between the wells. For a time- 




FIG. 2: Dependence of the final adiabatic population P of 
the 2nd well on the detuning rate a and nonlinearity A within 
the LZ (left) and GLZ (right) protocols. 



dependent coupling, the adiabatic following takes place 
under the condition [351 ] 

i|n(t)A(t) - n(t)A(t)| « \n 2 (t) + A 2 (t)\ 3 / 2 (32) 

which requires slow rates and large values of the coupling 
and/or detuning. 

Figure 1 shows that in the linear case both LZ and 
GLZ protocols give very similar results. The only tiny 
difference is that GLZ adiabatic energies fi±, calculated 
with (|23|) . less deviate from Ex,s at early and late time. In 
both protocols, the transfer is determined by the region 
where |/x_ — ~ fl(t — 0) = 1. As shown below, much 
more serious difference between LZ and GLZ appears in 
the nonlinear case. 



IV. RESULTS AND DISCUSSION 

The main results of our calculations are presented in 
Figs. 2-8. 

In Fig. [2j the final LZ and GLZ populations of the 
second well, P — A^(£ = +oo), are shown as a function 
of a for different values of the nonlinear parameter A. 
In all the cases, the BEC is initially in the first well, i.e. 
Ni(t = -oo) = l,N 2 (t = -oo) = 0. At the first glance, 
the LZ and GLZ transfers look quite similar. Without the 
interaction (A = 0) the processes are identical. For the 
repulsive interaction A = 4 and 20, both protocols pro- 
duce a complete transfer within a wide a-plateau. The 
stronger the interaction, the wider the plateau. At suffi- 
ciently high a, the process becomes too rapid and adia- 
batic transfer naturally fails, leading to decreasing P. In 
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Constant coupling 



Gaussian coupling 



1.0 

0.5 
N - 0.0 

r -t 

1.0 
0.5 
0.0 



A=4, a=3 



0.5 



40 



0.0 



-40 



A=-4,a=3 0.5 

0.0 



40 



-40 



t 



t 



FIG. 3: LZ (left) and GLZ (right) time evolution of popu- 
lations Ni (solid curve) and N2 (dash curve) of the first and 
second wells for the repulsive (a, b) and attractive (c, d) BEC. 




FIG. 4: GLZ adiabatic energies fi± (solid curves) for attrac- 
tive A=-4 (a) and repulsive A— 4 (b) interaction. The diabatic 
energies (dash lines) are given for the convenience of compar- 
ison. In both plots the rate a=3 is used. BEC is initially 
placed at the lower level as depicted by a bold dot. In the 
plot a) the terminal loop point and directions of tunneling 
to diabatic (arrow up) and adiabatic (arrow down) paths are 
marked. 



this region, behaviour of P depends on A but is the same 
for LZ and GLZ. 

It is remarkable that, unlike the STIRAP case [26[, the 
repulsive interaction (A > 0) not spoils but even favors 
the transport by forming the a-plateau. The robustness 
and completeness of the transport is illustrated in Fig. 
|3ja, b) for the particular case of A=4 and a=3. Instead, 
following Fig[^g-j) and Fig. [3Jc, d), the attractive in- 
teraction (A < 0) damages the adiabatic transport and 
the stronger the interaction, the less the final population 
P [45j |. So, in accordance with [l6|, we obtain the 
asymmetric interaction effect when the transfer crucially 
depends of the interaction sign. 

The origin of this effect is explained in Fig. [4] in terms 
of GLZ stationary eigenvalues \x± for the repulsive and 



2 5 Constant coupling 
a 



Gaussian coupling 
b) 




FIG. 5: LZ (left) and GLZ (right) chemical potentials fi± 
(solid curves) at different nonlinearity A and detuning rate 
a. The diabatic energies (dash lines) are given for the con- 
venience of comparison. In the plot c), the loop length limits 
(—tc.tc) are marked. 



attractive BEC. In both cases BEC is initially placed at 
the lower level, i.e. Ni(t — — oo) = l and A^(f = — oo)=0. 
The nonlinearity causes the loops in fi + (A=-4) and /i_ 
(A=4) levels. For A=-4, the loop takes place at the BEC 
transfer path and so a part of the condensate arrives to 
the loop terminal point. There is no way further and 
so BEC is enforced to tunnel to upper or lower levels as 
indicated by arrows in the plot a) . This results in leaking 
the population to the upper level /i_ and hence depletion 
of the adiabatic transfer. Instead, for A=4, the transfer 
path does not meet this obstacle and so is robust. 

The previous analysis mainly concerned the common 
features of LZ and GLZ protocols. Let's now consider 
the difference between these two cases, which takes place 
at small rates a. As seen from Fig. [2J for both protocols 
the nonlinearity results in a window < a < ot s where 
the adiabatic transfer completely vanishes. However, the 
GLZ window is much wider than the LZ one, which is a 
consequence of a finite duration of the Gaussian coupling. 

This point is scrutinized in Fig. [5] where behaviour of 
nonlinear loops in different transport regimes is demon- 
strated. In the first regime (plots (a,b)), the loops are 
absent because of a small nonlinearity A < 1. The LZ 
and GLZ stationary spectra and transport features look 
similar. In particular, in both cases the adiabatic trans- 
port is possible only in a narrow interval of small rates 
a [4q |. In the second regime with a larger nonlinearity 
A > Cl (plots (c-f)), the stationary spectra gain the loops 
but their length is still less than the Gaussian width, i.e. 
2t c <FWHM or t c < T. Then the LZ and GLZ spectra 



FIG. 6: GLZ evolution of adiabatic energies /i± (left) and 
population N2 (right) for the detuning rates a=0.328 and 
0.348. In all the panels A=4 is used. The vertical arrows 
indicate the minimal energy gap between /j,- (upper curve) 
and loop edge of (lower curve). 



u.u- 

1 2 3 4 5 

a 

FIG. 7: The adiabatic population P(a) for A =4 and Gaus- 

si 

ti 



remain to be similar and both protocols give a robust 
adiabatic transport P w 1 within common a intervals 
inside the plateau in Fig. 2. Note that according to (|29|) . 
the loop length 2t c rises with A and decreases with a as 
illustrated in Fig. 5. In other words, the loop behavior 
is determined not only by the nonlinearity but the de- 
tuning rate as well. The latter becomes decisive in the 
third regime that determines the P ~ windows at low 
rates a in Fig. 2. In this regime (plots (g-j)), the loop 
becomes comparable or longer than the coupling width, 
i.e. t c > r. Then the edges of the Gaussian coupling 
come to play and a new regime, which is absent in LZ, 
arises. At the edges we have f2(t) -C 1. So, in accor- 
dance to (|29|). the GLZ loop is extended and becomes 
longer than the LZ one, see plots (g, h). What is im- 
portant for our aims, in this regime the gap between /z_ 
and //+ shrinks at the loop edges, thus allowing a partial 
diabatic population transfer from /i + to /i_ and corre- 
sponding break of the adiabatic following. So, unlike the 
LZ, the main effect here takes place not at the coupling 
maximum but at its edges. For even longer loops (plot 
j)), the GLZ leads to merging and /z+ levels at early 
and late times, which fully breaks the adiabatic following 
and causes the windows exhibited in Fig. 2. Conversely, 
the constant coupling prevents the merging and so keeps 
the adiabatic transfer. 

The previous analysis suggests that at t c > T the suc- 
cess of adiabatic following is determined by the minimal 
gap between /i + and /Lt_. Fig. [6] illustrates this point 
for two small detuning rates which hinder and favor the 
transfer by giving iV2 ~0.02 and ~ 1, respectively. This 
example shows that Gaussian coupling provides a com- 
plete transfer switch by a tiny tuning of a and his switch 
is determined by the lateral gaps between /i + and 
This property can be used for effective coherent control 
of BEC transport. 



The Gaussian coupling delivers a new control param- 
eter T. Dependence of the transport on T is illustrated 



a) A > a > 



b) A < a < 



Il> |2> 



c) A > a < 




d) A < a > 




FIG. 8: Transport schemes for repulsive (A > 0, black dot) 
and attractive (A < 0, grey dot) BEC in directions indicated 
by arrows. Note relation of the sign of a with initial relative 
positions of the well depths. 



in Fig. \7\ It is seen that increasing T makes the Gaus- 
sian coupling closer to LZ case and leads to the downshift 
of the left plateau edge. Instead, decreasing r upshifts 
the plateau edge. Hence by changing T one can shift 
the switch value a s which, following a simple estimation, 
reads 



A 

2T ' 



(33) 



The figure shows that already at F < 3 the coupling 
becomes too short to support a slow adiabatic evolution 
and transport tends to vanish. Note that the value of T 
influences only the left side of the plateau and does not 
affect its right side. 

In the discussion above, the successful left-to-right 
transport of repulsive BEC (A > 0) was considered. This 
protocol is depicted in FigJS^,). However, by using the 
symmetry (|2ip . the similar transport can be produced 
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for the attractive condensate (A < 0) if we change the 
sign of the detuning rate a (FigJSJa)). Further, follow- 
ing the symmetry (|22|) . the right-to- left transport is also 
possible for the attractive (at a > 0) and repulsive (at 
a < 0) BEC, as is shown in Fig[8j;-d). So our protocol 
is quite universal. It can be also used for other kinds of 
multicomponent BEC, e.g. for two-component BEC in a 
single- well trap. 

As was mentioned in the previous sections, the GLZ 
is also a generalization of the Rosen-Zener method [S(| 
where the coupling is time-dependent but diabatic ener- 
gies are constant. Note that, unlike the GLZ, the non- 
linear RZ protocol totally blocks the adiabatic transfer 

M- 

Finally note that a possible (though weak) diabatic 
leaking Plz, appearance of the loops in nonlinear sta- 
tionary spectra, and finite rates a signify that the ac- 
tual transport is not perfectly adiabatic. However, be- 
ing driven by adiabatic protocols, the transport demon- 
strates high effectiveness and completeness in a wide 
range of nonlinearity. This proves that adiabatic popula- 
tion transfer schemes are indeed very robust and promis- 
ing. 

V. CONCLUSIONS 

We propose a simple and effective adiabatic population 
transfer protocol for the complete and irreversible trans- 
port of Bosc-Einstein condensate (BEC) in a double-well 
trap by using a pulsed coupling of a Gaussian form. Being 
mainly based on the familiar Landau-Zener (LZ) scheme 



with a constant coupling and having much in common 
with that scheme, our protocol is, nevertheless, more re- 
alistic and flexible. The protocol delivers the additional 
control parameter, the coupling width T, which allows to 
switch adiabatic transport at the detuning rate a beyond 
the LZ adiabatic limit a — > 0, and, what is important, 
provides a new transport regime when not the center but 
edges of the coupling are decisive. As a result, an effec- 
tive control of BEC transport, e.g. its rapid complete 
switch, becomes possible. In spite of asymmetric impact 
of nonlinearity, the protocol is quite universal and can 
be applied to both direct and inverse transport of BEC 
with both repulsive and attractive interaction. In other 
words, the protocol can always be cast to enforce the 
nonlinearity not hamper but favor the adiabatic trans- 
fer. The stronger the nonlinearity, the wider the range of 
the detuning rate a where the complete transfer is pos- 
sible. The protocol is actually a a generalization of both 
Landau-Zener and Rosen-Zener schemes. 

Our findings can open new opportunities in BEC dy- 
namics like excitation of topological states [HI, [3] , ex- 
ploration of diverse geometric phases generated in BEC 
transport jUj, etc. The later is especially interesting 
in relation to perspectives of using geometric phases in 
quantum computing [32], [H, [Hj] . 
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